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Abstract 

The deformation problem for a transversely isotropic elastic layer bonded to a rigid 
substrate and coated with a very thin elastic layer made of another transversely 
isotropic material is considered. The leading-order asymptotic models (for com¬ 
pressible and incompressible layers) are constructed based on the simplifying as¬ 
sumptions that the generalized plane stress conditions apply to the coating layer, 
and the flexural stiffness of the coating layer is negligible compared to its tensile 
stiffness. 
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1 Introduction 


Some natural biological tissues such as articular cartilage possess an inhomogeneous, lay¬ 
ered structure with anisotropic material properties. In particular, morphological studies of 
adult articular cartilage pHU] show three different zones of preferred collagen fiber bundle 
orientation. The superficial zone formed by tangentially oriented collagen fibrils provides 
a thin layer with a high tensile stiffness in the direction parallel to articular surface. It 
was shown Hi that the high transverse stiffness of the superficial tissue layer (character¬ 
ized by tangentially oriented collagen fibrils) is important in controlling the deformation 
response of articular cartilage. Generally speaking, the surface layer in biomaterials usu¬ 
ally has different mechanical properties than the underlying bulk material. Due to this 
circumstance, the mechanical deformation behavior is strongly influenced by the complex 
interaction of these layers [IB] . 

In the present paper, we consider the deformation problem for a transversely isotropic 
elastic layer reinforced with a thin elastic membrane ideally attached to one surface, 
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while the other surface is bonded to a rigid substrate. Following PEZ], it is assumed that 
the reinforcing layer is very thin (with respect to a characteristic size of the applied load) 
so that its deformation can be treated in the framework of the generalized plane stress 
state. Moreover, it is assumed that the flexural stiffness of the coating layer is negligible 
compared to its tensile stiffness. Thus, the reinforcing thin layer is regarded as an elastic 
membrane. 

The rest of the paper is organized as follows. In Section [2j we formulate the three- 
dimensional boundary conditions for a coated elastic layer. The deformation problem 
formulation itself is given in Section [3j Asymptotic analysis of the deformation problem is 
presented in Section |4} The main result of the present paper is presented by the leading- 
order asymptotic models for the local indentation of the coated clastic layer developed 
in Section [5] for the cases of compressible and incompressible layer. Finally, Section [6] 
contains some discussion of the obtained asymptotic models and outlines our conclusion. 


2 Boundary conditions for a coated elastic layer 


We consider a very thin transversely isotropic clastic coating layer (of uniform thickness h) 
ideally attached to an clastic layer (of thickness h) made of another transversely isotropic 
material (see Fig. [T|. Let the five independent elastic constants of the elastic layer and its 
coating are denoted by An, A 12 , A 43 , A 33 , A 44 and An, A 12 , Ai 3 , A 33 , A 44 , respectively. 


h 0 
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Figure 1. An elastic coated layer as a model for articular cartilage (the histological image of 
articular cartilage is taken from the paper EE]). The regions —h<z< 0 and z > h represent 
the superficial zone and the subchondral bone, respectively. 

Under the assumption that the two layers are in perfect contact with one another along 
their common interface, z = 0, the following boundary conditions of continuity (interface 
conditions of perfect bonding) should be satisfied: 

v(y,0) = v(y,0), w>(y, 0) = w(y, 0), (1) 

&3j(y, o) = 03 j(y, o), j = i, 2,3. (2) 

Here, (v, w) is the displacement vector of the elastic coating layer, a tJ are the correspond¬ 
ing components of stress. In what follows, we make use of the Cartesian coordinate system 
(y ,z), where y = ( 2/1 , 2 / 2 ) are the in-plane coordinates. 

On the upper surface of the two-layer system, z = —h, we impose the boundary conditions 
of normal loading with no tangential tractions 


0-31 (y, ~h) = cr 32 (y, -h) = 0, <j 33 (y, -h) = -p( y), 


( 3 ) 










where p(y) is a specified function of external loading. 

Following Rahman and Newaz [IS], we simplify the deformation analysis of the elastic 
coating layer based on the following two assumptions: (1) the coating layer is assumed to 
be very thin, so that the generalized plane stress conditions apply; (2) the flexural stiffness 
of the coating layer in the z-direction is negligible compared to its tensile stiffness. 


In the absence of body forces, the equilibrium equations for an infinitesimal element of 
the coating layer are 


doji , d°j3 

dyi dy 2 dz 


0, j = 1,2,3. 


(4) 


The stress-strain relationship for the transversely isotropic elastic coating layer is given 
by 
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where 2A 66 = An — A 12 . 


Integrating Eqs. (J4]) through the thickness of coating layer and taking into account the 
interface and boundary conditions ([2]) and (|3]) , we get 
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Here, o l3 are the averaged stresses, i.e., 


( 6 ) 

(7) 



y,z) dz. 


Under the simplifying assumptions made above, we have 

^13 = ^23 = 033 = 0- (8) 

Hence, Eq. ([7]) immediately implies that 
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(9) 


















Moreover, in view of (j8|, the averaged strain i 33 must satisfy the equation 

A\ 3 £\\ + A 13 £ 2 2 + ^433^33 — 0. 


Therefore, the in-plane averaged stress-strain relationship takes the form 
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where we introduced the notation 
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On the other hand, in view of the interface conditions ([Tj) , we will have 
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where £n, £ 22 , and £12 are the in-plane strains in the coated elastic layer 0 E (0, h ). 


Therefore, taking Eqs. (10) and (12) into account, we transform the boundary conditions 
@ as follows: 
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Finally, the above boundary conditions can be rewritten in the matrix form as 


C31 e l + 032^2 
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where £(V y ) is a 2 x 2 matrix differential operator such that 
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Thus, the deformation problem for an elastic layer coated with a very thin flexible elastic 
layer is reduced to that for the elastic layer without coating, but subjected to a different 
set of boundary conditions ([9]) and (13) on the surface z = 0. 





















3 Deformation problem formulation 


Now, let us consider a relatively thin transversely isotropic elastic layer of uniform thick¬ 
ness, h, coated with an infinitesimally thin clastic membrane and bonded to a rigid sub¬ 
strate, so that 
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(15) 


In the absence of body forces, the vector (v, w) of displacements in the elastic layer satisfies 
the Lame system 


<9 2 v d 

+ ( A u — A m )W yS7 y - v + A 44 + (^13 + ^ 44 ) y w = 

d 2 w d 

A 44 Aj / u' + A 33 -^-^- + (A i3 + A44)— V y - v = 0. 


(16) 


Assuming that the coated layer is loaded by a normal load and denoting the load density 
by p, we require that 

= -p. (17) 
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Based on the analysis performed in Section [2j the influence of the elastic membrane 
(coating layer) is introduced by the boundary condition 
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where £(V y ) is the matrix differential operator dehned by formulas (14) 
Taking into account the stress-strain relationship 
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we rewrite Eqs. (17), (18) as follows: 
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( 20 ) 


Equations (15), (16), (19), and (20) comprise the deformation problem for the coated 


transversely isotropic elastic layer bonded to a rigid substrate. 

























4 Asymptotic analysis of the deformation problem 


Let /i* be a characteristic length of the external load distribution. Denoting by e a small 
positive parameter, we require that 

h = eh * (21) 

and introduce the so-called “stretched” dimensionless normal coordinate 

C = —. 
eh* 

In addition, we non-dimensionalize the in-plane coordinates by the formulas 

Vi=^r, * = 1 , 2 , rj — (jji, 772 ), 

/ 1 * 


so that 


<9 _ i a _ 1 

<9z eKdC’ y h* 


Moreover, we assume that the tensile stiffness of the coating layer is relatively high, i.e., 
A u >>> A 11; and so on. To fix our ideas, we consider the situation when 


£(v„) = r‘r(v,), 


( 22 ) 


that is, in particular, the ratio An/An is of the order of e. 

Employing the perturbation algorithm |9j, we represent the solution to the deformation 


problem (15), (16), (19), (20) as follows: 


v = £ 2 v 1 (r/,C) + ..., 

W = £W°(r), C) + e 3 W 2 (r], 0 + ... . 

Here, for the sake of brevity, we include only non-vanishing terms. 
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It can be shown (see, in particular, 0) that the leading-order term in rt24j) is given by 

w°(v,0 = h n(l-(), (25) 

2^33 


whereas the first non-trivial term of the expansion (23) satisfies the problem 
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Substituting the expansion (25) for w° into the above equations, we obtain 
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The solution to the boundary-value problem (26) is represented in the form 


2 ^ 33^44 

where V 1 (r/) satisfies the equation 
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For the second non-trivial term of the expansion (24), we derive the problem 
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(28) 


Now, substituting the expressions (25) and (27) for w° and v 1 , respectively, into the above 
equations, we arrive at the problem 
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Integrating Eq. (29) twice with respect to ( and taking into account the boundary condi¬ 
tion (30)2, we obtain 
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where C' 2 (t 7 ) is an arbitrary function. 




































The substitution of (31) into the boundary condition (30h yields 
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Hence, in view of this relation, formula (31) implies 
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where V 1 is the solution of Eq. (28) 


5 Local indentation of the coated elastic layer: Leading-order asymptotics 
for the compressible and incompressible cases 


Recall that the local indentation of an elastic layer is defined as 


w 0 (y) = w(y>°)> 

where w(y, 0 ) is the normal displacement of the layer surface. 


In the case of compressible layer, Eqs. (24) and (25) yield 


^o(y) - u—p(y), 
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(33) 


so that the deformation response of the coated elastic layer is analogous to that of a 
Winkler elastic foundation [T31IT6] with the foundation modulus k = A 33 // 1 . In other 
words, the deformation of the elastic coating does not contribute substantially to the 
deformation of a thin compressible layer. 


When the layer material approaches the incompressible limit, the right-hand side of (33) 
decreases to zero and the first term in the asymptotic expansion (24) disappears. At that, 
based on the known results HU, it can be shown that the ratios A 13 /A 33 and A 44 /A 33 
tend to 1 and 0 , respectively. 


Therefore, in the incompressible limit situation formula (32) reduces to 
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where a 44 = A 44 is the out-of-plane shear modulus of the elastic layer, while V 1 (r/) satisfies 
the equation 

1 ^ . h 
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Thus, in the case of incompressible bonded elastic layer, formulas (22)—(24), (34), and 
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where the vector v 0 (y) solves the equation 
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(37) 


Observe that in view of (27), the vector-function v 0 (y) has the meaning of the tangential 


displacement of the surface point (y,0) of the clastic layer. Recall also that the matrix 
differential operator £(V y ) is defined by formulas (14). 


6 Discussion and conclusion 


Observe that in the axisymmetric case (in view of ( 12 )), we have 

a rr = Au£ rr + Ai2Eoei o'ee — A^rr + Ai\£qo, o r e = 0 , 

where 
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while Eqs. (J 6 ]) should be replaced with the following one: 
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Correspondingly, the boundary condition (13) takes the following form: 
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ft is to note here that the axisymmetric boundary condition (38) was previously derived 
in a number of papers [TlGiTTITS] . 

Finally, let us consider two opposite limit situations in the incompressible case. First, 
when the coating is absent and £(V :y ) = 0, Eq. (37) implies 


v 0 (y) = 
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The substitution of this expression into formula (36) gives 
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which completely agrees with the asymptotic model developed in (see also mra. 

Second, in the case of very stiff (inextensible) coating we will have v 0 (y) = 0, and formula 
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(36) reduces to 
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In other words, comparing (39) and (40), we conclude that the inextensible membrane 
coating attached to the surface of a thin incompressible elastic layer reduces the out-of¬ 
plane shear compliance of the bonded layer by four times. 


The main result of the paper, presented in Section [4j is the leading-order asymptotic 
formulas for the normal displacement (called the local indentation) of the surface points 
for compressible and incompressible thin bonded elastic layers reinforced with an elastic 
membrane. 


7 Acknowledgment 


The authors are grateful for support from the FP7 IRSES Marie Curie grant TAMER 
No 610547. 


References 


[1] V.M. Alexandrov, S.M. Mkhitaryan, Contact Problems for Solids with Thin Coatings and 
Layers [in Russian], Nauka, Moscow, 1985. 

[2] 1.1. Argatov, Pressure of a paraboloidal die on a thin elastic layer, Doklady Phys. 50 (2005) 
524-528. 

[3] 1.1. Argatov, Development of an asymptotic modeling methodology for tibio-femoral contact 
in multibody dynamic simulations of the human knee joint, Multibody Syst. Dyn. 28 (2012) 
3-20. 

[4] I. Argatov, G. Mishuris, Elliptical contact of thin biphasic cartilage layers: Exact solution 
for monotonic loading, J. Biomech. 44 (2011) 759-761. 

[5] I. Argatov, G. Mishuris, Frictionless elliptical contact of thin viscoelastic layers bonded to 
rigid substrates, Appl. Math. Model. 35 (2011) 3201-3212. 

[6] V.I. Avilkin, V.M. Alexandrov, E.V. Kovalenko, On using the more-accurate equations of 
thin coatings in the theory of axisymmetric contact problems for composite foundations, J. 
Appl. Math. Mech. 49 (1985) 770-777. 

[7] J.R. Barber, Contact problems for the thin elastic layer. Int. J. Mech. Sci. 32 (1990) 129-132. 

[8] I.C. Clarke, Articular cartilage: a review and scanning electron microscope study. 1. The 
interterritorial fibrillar architecture, J. Bone Joint Surg. Br. 53-B (1971) 732-750. 

[9] A.L. Gol’denveizer, Derivation of an approximate theory of bending of a plate by the method 
of asymptotic integration of the equations of the theory of elasticity, J. Appl. Math. Mech. 
26 (1962) 1000-1025. 

[10] L.C. Hughes, C.W. Archer, I. ap Gwynn, The ultrastructure of mouse articular cartilage: 
collagen orientation and implications for tissue functionality. A polarised light and scanning 
electron microscope study and review, Europ. Cells Materials 9 (2005) 68-84. 




[11] M. Itskov, N. Aksel, Elastic constants and their admissible values for incompressible and 
slightly compressible anisotropic materials. Acta Mech. 157 (2002) 81-96. 

[12] M.J. Jaffar, Asymptotic behaviour of thin elastic layers bonded and unbonded to a rigid 
foundation. Int. J. Mech. Sci. 31 (1989) 229-235. 

[13] K.L. Johnson, Contact Mechanics. Cambridge Univ. Press, Cambridge, UK, 1985. 

[14] R.K. Korhonen, M. Wong, J. Arokoski, R. Lindgren, H.J. Helminen, E.B. Hunziker, J.S. 
Jurvelin, Importance of the superficial tissue layer for the indentation stiffness of articular 
cartilage, Med. Eng. Phys. 24 (2002) 99-108. 

[15] R.V. Patel, J.J. Mao, Microstructural and elastic properties of the extracellular matrices of 
the superficial zone of neonatal articular cartilage by atomic force microscopy, Front. Biosci. 
8 (2003) al8-25. 

[16] V.L. Popov, M. HeB, Methode Der Dimensionsreduktion in Der Kontaktmechanik Und 
Reibung [in German], Springer, Berlin, 2013. 

[17] M. Rahman, G. Newaz, Elastostatic surface displacements of a half-space reinforced by a 
thin film due to an axial ring load, Int. J. Eng. Sci. 35 (1997) 603-611. 

[18] M. Rahman, G. Newaz, Boussinesq type solution for a transversely isotropic half-space 
coated with a thin film, Int. J. Eng. Sci. 38 (2000) 807-822. 



